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SUMMARY 


A  new  time  series  model  for  exponential  variables  having  first 
order  autoregressive  structure  is  presented.   Unlike  the  recently  studied 
standard  autoregressive  model  in  exponential  variables  (EAR(l)),  runs  of  con- 
stantly scaled  values  are  avoidable,  and  the  two  parameter  structure  allows  some 
adjustment  of  time  nonreversibility  effects  in  sample  path  behavior.   The 
model  is  further  developed  by  the  use  of  cross-coupling  and  antithetic 
ideas  to  allow  negative  dependency.   Joint  distributions  and  autocorrelations 
are  investigated.   A  transformed  version  of  the  model  has  a  uniform  marginal 
distribution  and  its  correlation  and  regression  structures  are  also  obtained. 
Estimation  aspects  of  the  models  are  briefly  considered. 


KEYWORDS:   Autoregressive  model  in  exponential  variables;  Negative  correlation; 
Cross-coupled  processes;  Antithetic  variables;  Correlated  uniform 
process:  Time  series;  Point  process;  Simulation. 


1 .   INTRODUCTION 

In  this  paper  we  begin  by  introducting  a  new  two-parameter  model, 
to  be  called  NEAR(l),  first  mentioned  in  Lawrance  (1979),  for  a  first-order 
autoregressive  time  series  with  exponentially  distributed  marginals.   The 
model  is  a  first-order  Markov  process.   Suitably  choosing  one  of  the 
parameters  as  a  function  of  the  other  produces  a  one-parameter  first- 
order  autoregressive  process  which  can  give  any  value  of  the  lag  one 
autocorrelation  between  zero  and  one.   One  particular  model  produced  in 
this  way  is  the  EAR(l)  model  introduced  by  Gaver  and  Lewis  (1980) ;  this 
model  had  the  problem  that  a  "zero-defect"  caused  successive  values  of 
the  process  to  be,  at  times,  fixed  multiples  of  the  previous  values.   The 
NEAR(l)  model  does  not  have  this  defect  except  for  the  EAR(l)  special  case 
and  thus  seems  much  more  suitable  than  the  EAR(l)  model  for  the  modelling 
of  real  data.   In  addition,  the  fact  that  there  are  two  parameters 
indexing  the  dependency  structure  of  the  model  allows  one  to  consider 
sample  path  behavior  as  well  as  the  customary  fitting  of  the  first  and 
second  order  moments  to  the  data.   The  model  is  defined  in  Section  2. 

At  another  extreme  from  the  EAR(l)  model,  a  one-parameter  model 
(TEAR(l))  is  produced  which  is  much  easier  to  extend  to  higher  order 
autoregressive  structures  than  is  the  EAR(l)  model  (Lawrance  and  Lewis, 
1980).   However  while  it  has  no  zero  defect,  this  TEAR(l)  model  produces 
realizations  which,  for  high  serial  correlation,  tend  to  run  up  most  of 
the  time;  for  the  general  NEAR(l)  model  these  aspects  can  be  adjusted. 
A  one-parameter  model  which  can  mimic  some  of  the  time-reversible  char- 
acter of  normal  AR(1)  processes  is  produced  from  the  NEAR(l)  model  by 
requiring  either  that  the  probability  of  a  jump  up  from  one  value  to 
the  next  be  one-half  or  requiring  that  the  first  directional  moments  be 
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equal.   A  property  which  the  NEAR(l)  model  does  not  share  with  its  special 
EAR(l)  case  is  additivity,  so  that  extensions  to  Gamma  marginals  are  not 
automatic;  other  marginal  distributions  are  possible  with  the  NEAR(l) 
structure  but  these  are  not  discussed  here. 

An  important  property  of  the  NEAR(l)  models  is  that  they  are  simple 
random  linear  combinations  of  independent  exponential  variables  and  there- 
fore easy  to  simulate.   This  simplicity  is  bought  at  the  price  of  auto- 
correlations which  are  nonnegative. 

The  second  thrust  of  the  paper  concerns  alternation  and  negativity 
of  autocorrelations;  this  will  be  achieved  by  a  scheme  coupling  two  antithetic 
NEAR(l)  sequences,  a  scheme  introduced  by  Gaver  and  Lewis  (1980)  for  the  negatively 
correlated  EAR(l)  process.   The  resulting  model,  to  be  called  the  NEARA(l) , 
includes  both  the  NEAR(l)  and  hence  TEAR(l)  as  special  cases;  it  has  auto- 
correlations which  alternate  into  negativity  under  a  geometrically  decaying 
envelope.   However,  simulation  of  the  negatively  dependent  models  involves 
random  linear  combinations  from  independent  pairs  of  negatively  dependent 
exponential  variables,  and  this  can  be  complicated.   Most  developments  in 
the  paper  are  undertaken  for  the  general  NEARA(l)  model,  and  further  detail- 
ing of  results  are  given  separately  for  the  positive  and  negative  dependency 
cases.   In  particular,  the  paper  deals  with  the  allowable  range  of  lag  one 
autocorrelations,  lag  r  bivariate  distributions,  exponentiation  of  the 
models  to  have  uniform  marginal  distribution,  and  aspects  of  time  reversi- 
bility, sample  path  behavior  and  estimation. 

Simulation  aspects  of  the  models  are  discussed  in  Lawrance  and  Lewis 
(1980);  detailed  graphical  representations  of  different  sample  path 
behaviors  are  also  given  there. 


2.   CONSTRUCTION  OF  THE  MODELS 

The  conventional  linear  autoregressive  model  (AR(1))  with  exponential (X) 
marginal  distributions  (Gaver  and  Lewis,  1980)  takes  the  form 


X  =  pX    +  i 
n     n-1 


0       w.p .   p 

n  =  0,1,2,...  ,       (2.1) 
E       w.p.  1-p 


where  p   is  a  parameter   (0  ;!  p  <  1)   and  the   E  ,  n  =  0,1,2,...   are  inde- 
pendent exponential  variables  with  parameter  X  >  0.   This  EAR(l)  model  has 

serial  correlations  of  order  r,  p   =  corr(X  ,X  ,  ) ,  given  by  p    and  gen- 

r         n  n+r   °      J  ° 

erates  sample  paths  in  which  large  values  are  followed  by  runs  of  falling 

values  with  geometrically  distributed  run-length.   The  large  values  arise 

when  E   is  included,  while  the  falling  values  stem  from  the  selection  in 
n 

(2.1)  giving  only  X  =  pX  . .   This  behavior  is  likely  to  limit  the  broad 

n     n-1 

applicability  of  the  model,  although  it  can  be  overcome  the  more  complicated 
moving-average  and  mixed  moving  average-autoregressive  developments  (Lawrance 

and  Lewis,  1977,  1980a;  Jacobs  and  Lewis,  1977). 

An  alternative  exponential  first-order  autoregressive  Markov 
model  is  obtained  by  interchanging  the  independent  and  identically  distri- 
buted variables   X  n   and  E   in  (2.1);  this  can  have  no  effect  on  the 

n- 1  n 

exponential(X)    marginal  distribution  of     X    's.      Proceeding   this  way,   with 
p      replaced  by      1-ot ,    we  have    the   model 


X      =    (l-a)E   +< 
n  n 


X     .,        w.p .   a 
n-1 

n   =   0.k,2, ...    .  (2.2) 

0  w.p.    1-a 


This  exponential  AR(1)  model,  called  TEAR(l) ,  is  again  Markovian  and  has 
the  a   correlation  structure  of  the  EAR(l)  model;  it  is,  as  will  be  shown 


later,  particularly  tractable  analytically.   The  characteristic  behavior  of 

realizations  generated  by  this  model  (particularly  distinct  when  a  is 

large)  is  that  of  runs  of  rising  values  (with  geometrically  distributed 

run  length)  when  the  selection  (l-a)E  +  X  1   is  being  made,  followed  by 
0  n    n-1 

a  sharp  fall  when  the  selection   (l-a)E    is  made  without  inclusion  of  the 

n 

previous  value.   Illustrations  of  these  effects  both  for  the  EAR(l)  and 

TEAR(l)  models  are  given  in  the  simulations  of  Fig.  la  and  Fig.  lb.   These 

simulated  sample  paths  use  the  same  simulated  exponential  error  sequence 

(E  )  . 
n 

Broader  behavior  in  realizations  generated  by  an  exponential  model 

can  be  obtained  from  the  model  in  which  the   X   .,   of  (2.2)   is  scaled  by 

n-1  J 

a  coefficient  6.   This  gives  the  proposed  NEAR(l)  model  (Lawrance ,  1980) 


as 

I 


X  =  e   +  < 
n    n 


BX  .    w .p.  a 
n-1 


n  =  0,1,2,...  ,         (2.3) 
w.p.   1-a 


where  the  existence  and  distribution  of  the  i.i.d.  {e  }   sequence  which 

makes  the   X  ' s  in  the  stationary  case  have  exponential (X)  distributions, 

needs  to  be  established  afresh.   We  now  show  that   e   must  have  a 

n 

particular  mixed  exponential   distribution. 

Let   the   Laplace-Stieltjes   transforms   of   the      X     and     e       variables 
be   denoted  by 

<|>x(s)    =   E(e"   X}  and  <J.    (s)    =   E{e"Se}    .  (2.4) 


Then  (2.3)  gives,  if  we  assume  stationarity , 


i  t     \     X A  T  PS  A ,  «  j-n 

Vs;  !  a<b  (Bs)  +  (1-a)  :  X  +  s   X  +  (l-a)Bs   '  K        J 

x 


on  using  <J>(s)  =  X/(X  +  s).   Thus,  providing  a   and   B   are  not  both  equal 
X 

to  one,   e   can  be  generated  from  an  E   by  the  exponential  mixture 
n  n   j  *■ 


e   = 
n 


T7  !  -  & 

E  w.p 


n  ,F-   1  -  (l-a)B 

n  =  0,1,  ...  .      (2.6) 

(1"a)BEn      "•"•   1  -  (l-a)B 

When  a  =  0  or  B  =  0   the   (X  }   are  exponential  i.i.d.,  whereas  with 

n 

a  =  1   the  EAR(l)  model  (2.1)  is  obtained  with  p   =  B.   When  B  =  1   the 

TEAR(l)  model  is  obtained.   Thus  the  two-parameter  exponential,  first-order, 

autoregressive  Markov  NEAR(l)  model  can  be  expected,  for  fixed  serial 

correlation  of  lag  1,   p   =  aB,  to  model  broader  behavior  than  is  obtained 

in  the  extreme  cases  (a  =  1   or   B  =  1) .   In  particular  a  and  B   can  be 

chosen  to  produce  both  runs  of  ascending  and  descending  values,  intermediate 

to  the  profiles  of  EAR(l)  and  TEAR(l)  models,  as  was  illustrated  in  Fig.  la 

or  Fig.  lb.   Figure  lc  represents  an  intermediate  case  which  will  be  discussed 

in  Section  8.   Note  that  the  correlation  is  the  same,  0.75,  in  all  three  figures 

It  is  also  clear  irom  the  Markovian  nature  of  the  model  (i.e.  that 

conditional  on  X   n  =  x  n   the  distribution  of  subsequent  values 

n-1    n-1 

X  ,  X  .,  ...    is  independent  of  X   ~,  X   _.,...),  that  if  X_   is 
n   n+1  r  n-2   n-3  0 

exponential(X)  and  independent  of   E  ,  E  ,  ...  ,  then  the  process   X  , 

n  =  1,2,...   is  stationary.   Note  too  that  the  NEAR(l)  model  is,  by 

definition,  explicitly  (physically)  autoregressive  and  thus  not  only  auto- 

regeressive  in  the  sense  that   E(X  |x   ..  =  x)   is  a  linear  function  of   x. 

n   n-1 


We  note  too  that  the  NEAR(l)  model  gives  a  solution  to  the  (random) 
stochastic  difference  equation 


X  =  A  X 


n 


n  n-1 


+  B 


n 


n  =  0,1,2,.. 


(2.7) 


discussed  by  Vervaat  (1979)  in  which  k^  =   B  w.p.  a  and  A^  =   0  w.p.  (1-a): 
Vervaat's  paper  discusses  questions  of  existence  and  infinite  divisibility 
applying  to  the  model  (2.7). 

In  the  NEAR(l)  model  the  parameters  a   and  B   are  nonnegative. 

Therefore  the  autocorrelations   p   =  (aB)    are  positive  and  geometrically 

decreasing.   This  is  unlike  the  standard  AR(1)  model  with,  say,  normal 

marginals,  where   p    can  be  negative,  so  that  the  autocorrelations  can 

alternate  between  positive  and  negative  values  with  a  geometrically 

decreasing  envelope.   To  extend  the  exponential  models  to  the  situation  where 

there  is  a  possibility  of  alternation  in  the  autocorrelations  and  negative 

correlation  requires  some  sacrifice  of  simplicity.   As  noted  in  Section  1, 

the  primary  idea  here  is  to  cross -couple  two  sequences   (X  }   and   (X') 

n  n 

with  identically  exponentially  distributed  marginal  distributions  across 
an  independent  bivariate  sequence    {c  ,  e'}   of  negatively  correlated 
and  marginally  identical  variables.   This  final  development  produces  our 
so-called  NEARA(l)  model;  it  is  specified  by  the  equations 


X  =  e   +  BV  X»    , 
n    n     n  n-1 


X'  =  e*  +  BV'X  ,  , 
n    n     n  n-1 


V  = 

n 


1   w.p.  a 


0   w.p.  1-a, 


n  =  0,1,2, . ..   (2.8) 


1 

w.p . 

a 

V    = 
n 

0 

w.p. 

1-a 
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where  the  serially  independent  binary  pairs   V   and  V'   generally  have 

n        n 

negative  dependency.   Some  insight  into  the  model  comes  from  seeing  that 

X   is  positively  dependent  on  X'_  ;  this  is  negatively  dependent  on  X  , , 

so  making  X   and  X     negatively  dependent.   Though  defined  compactly 

in  terms  of  the  two  processes,  the  interest  here  is  in  the  marginal  process 

X  .   A  univariate  description  of  X   is  possible  and  given  at  equation 

(3.2). 

The  special  case  of  the  bivariate  sequences   {e  ,  e'}   and   {V  ,  V'} 

n   n  n   n 

in  which  e   =  e1   and   V  =  V'   recovers  the  NEAR(l)  model.   The  special 
n     n        n    n 

case  when  8=1  will  be  called  the  TEARA(l)  model. 

Detailed  aspects  of  the  sequence   {X  }   depend  on  the  joint  distri- 
butions of   {e  ,  e'}   and   {V  ,  V'},  though  the  marginal  distributions  of 
n   n  n   n 

e    and  e1   must  be  as  at  (2.6)  for  X   to  be  marginally  exponential.   For 
n        n  n  r 

instance  for  the  TEARA(l)  model,  strongest  alternation  in  serial  correlations 
is  obtained  when  the  {e  ,  e'}   are  maximally  negatively  correlated  expo- 
nential variables  and  therefore  are  antithetic  pairs,  and  similarly  for 
the  binary  pairs   {V  ,  V'}.   For  the  broader  NEARA(l)  model,  (2.8), 
negatively  correlated  mixed  exponential  variables  {e    ,  e'}   are  required. 
Some  of  these  aspects  of  the  model  are  explored  in  general  and  for  specific 

{e  ,  e1}   and   {V  ,  V'}   distributions  in  Sections  4  and  5.   In  this 
n   n  n   n 

respect  this  paper  extends  results  and  details  for  the  negatively  correlated 

EAR(l)  model  given  by  Gaver  and  Lewis  (1980) . 

Note  that  while   {X  ,  X'}   is  a  bivariate  Markovian  model,  the 

n   n 

full  Markovian  property  of  X   individually  is  lost  unless  it  reduces  to 

n 

the  NEAR(l)  model;  that  is  to  be  expected  from  the  cross-dependency  built 
into  the  model. 


3.   AUTOCORRELATION  STRUCTURE  OF  THE  MODEL 

The  simple  autocorrelation  structure  of  the  NEARA(l)  model,  as 
given  at  (2. 8)  by 


X  =  e   +  BV  X'    ,     X'  =  e'  +  BV'X  ,  (3.1) 

n    n     n  n-1        n    n     n  n-1 


is  best  approached  by  recursively  expressing  the  dependency  of  X   on 

either  X'   ,  X  „,  X'  _,  X   .,...,   and  so  on.   Directly  from  (3.1)  it 
n-1   n- I        n-J   n-H 

can  be  noted  that  the  distribution  of   (X  ,X  .)   is  simply  expressed  in 

n  n-1  J        r 

terms  of  the  distribution  of   (X  1  ,X'   ):   in  fact,  this  latter  joint 

n-1  n-1 

distribution,  equivalently   (X  ,X'),  plays  a  central  role  in  the  process. 

However,  substituting  for  X'  ,   in  the  first  equation  (3.1)  from  the 

n-1 


second,  gives 


X  =  £   +  BV  e'   +  bV  nV  X   .  .  (3.2) 

n    n     n  n-1      n-1  n  n-2 


Hence  the  joint  distribution  of   (X  ,X   _)   does  not  need  to  be  expressed 

n  n-z 

in  terms  of   (X  ,X')   and  this  is  very  convenient.   Generally,  there  is  this 
n  n 

distinction  between  the  odd-r  and  even-r  cases  of   (X  ,X   ) .   This  is  shown 

n  n-r 

in  the  following  key  expressions  which  are  obtained  by  repeated  substitutions; 

X  =  e   +  BV  e'   +  B2V  ,V  e   -  +  •••  +  (Br_1V'   ,  „  •••  V'  n£    n) 
n    n     n  n-1      n-1  n  n-2  n-r+2      n-1  n-r+1 

+  (BrV  ^,    ••■  V  .V  X1   );       (r  odd)  (3.3) 

n-r+1      n-1  n  n-r 


X     =  £      +  BV  e'        +  B2V    .V  £      _   +  •••    +(Br   1V        ,.   •••    V      V  £f        ,) 

n  n  n  n-1  n-1  n  n-2  n-r+2  n-1  n  n-r+1 


+  (eV      j_.    •••    V      V  X        ).  (r      even)  (3.4) 

n-r+1  n-1  n  n-r 


8 


The  autocovariances  of   {X  }   follow  easily  from  (3.3)  and  (3.4).   On 


n 


noting  that  the  indicator  variables   {V.}   occur  independently  in  the 


products,  we  have 


and  hence 


E(V   ,_ 
n-r+1 


V'   V  )  =  a' 
n-1  n 


(3.5) 


(a6)   Var(X   ) 
n-r 


Cov(X  ,X    )  ■  < 
n  n-r 


(r  even) 


(a3)   Cov(X   ,X'   )    (r  odd)  . 
n-r  n-r 


(3.6) 


In  terms  of  correlations,  this  central  result  becomes 


Corr(X  ,X   )  =  < 
n  n-r 


(aS) 


(a8)   Corr(X  ,X') 
n  n 


(r   even) 


(r   odd) 


(3.7) 


Alternation  of  these  autocorrelations  under  a  geometric  envelope  is  evident; 

negativity  of  the  odd  lag  correlations  requires  the  negativity  of 

Corr(X  ,X').   For  the  simpler  NEAR(l)  model  in  which  X  =  X1   the  Markov 
n   n  n    n 

(aB)    correlation  structure  is  evident. 

For  the  NEARA(l) ,  an  investigation  of   Corr(X  ,X')   is  required. 

n     n 

To    this   end,    multiply   together   the   respective   sides   of    the   two   equations 
(3.1),    giving 


X  X'    =   e    e'    +  BV'e   X      .    +     6V  e'X'    ,    IM  V'X     ,X'    , 
n  n  n  n  n  n  n-1  n  n  n-1  n  n  n-1   n-1 


(3.8) 


and  take  expectations.   Let 


€  =  Cov(e  ,e') 
n  n 


and 


v   =  Cov(V  ,V) 
n  n 


(3.9) 


and  assume   stationarity .      Then   following   from   (3.8)    there   is    the   result 


Corr(X   ,X')    -   (€  +  B2v)/{1  -    (a2+v)g2}. 
n     n 


(3.10) 


The  important  conclusion  is  that  maximum  negativity  of   Corr(X  ,X')   is 

n  n 

obtained,  for  any  fixed  values  of  a   and  6,  for  maximum  negative  correla- 
tions within  the  pairs   (e  ,e')   and   (V  ,V).   The  proof  is  omitted. 

n  n  n  n 

Obtaining  this  maximum  negative  correlation  by  the  use  of  antithetic 
variables  is  developed  in  the  next  section. 


4.   ANTITHETIC  ASPECTS  OF  THE  MODEL. 

It  is  simplest  to  deal  first  with  the  binary   (V  ,V)   variables; 

n  n 

the  basic  antithetic  idea  is  to  relate  the  distribution  of  V   to  a 

n 

monotonic  transformation  of  a  uniform  variable   U  on   (0,1);  then  V1 


n 


is  the  same  transformation  of   1  -  U  which  also  has  a  (0,1)  uniform 


distribution.   The  variables   V   and  V'   are  then  maximally  negatively 

n        n  j        o  j 

correlated.   Thus,  we  define 


V  =  1   if   U   <  a 
n  n  — 


V  =  0   if   U   >  a 
n  n 


V'  =  1   if  1-U  <  a   or   U   >  1-a 
n  n  —        n  — 


V  =  0   if   1-U   >  a  or   U  <  1-a 
n  n  n 


The  resulting  joint  distribution  takes  one  of  two  forms,  as  given  below, 
depending  on  whether  a  <  1/2   or  a  >  1/2: 
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V      = 
n 

1 

0 

V 

n 

V'    =   1 
n 

0 

a 

a 

0 

a 

l-2a 

1-a 

V 
n 

a 

1-a 

1 

V      = 
n 

1 

0 

V 

n 

V    =   1 
n 

2a- 1 

1-a 

a 

0 

1-a 

0 

1-a 

V 
n 

a 

1-a 

1 

(4.1) 


(a  <  1/2) 


(a  >  1/2) 


2  2 

The  resulting  covariances  are   v  =  -a   for  0  _<  a  _<  1/2   and  v  =  (1-a) 

for   1/2  <_  a  <  1;  the  corresponding  correlations  are  thus  -a/(l-a)   if 
0  _<  a  _<  1/2   and  -(1-a) /a   if   1/2  <_  a  <  1 .   The  a  =  1   case  is  exceptional 
and  is  excluded  since  it  leads  to  the  negatively  correlated  EAR(l)  model 
treated  in  Gaver  and  Lewis  (1980) . 

Next  we  consider  how  to  obtain  a  bivariate  mixed  exponential  dis- 
tribution for   (e  ,e')   having  maximum  negative  dependency.   In  the  case 
of  positive  continuous  random  variables,  the  maximum  negative  correlation 

is  obtained  by  the  antithetic  pair  (Moran,  1967).   However,  with   (e  ,e') 

n  n 


having  mixed  exponential  marginals,  the  full  antithetic  distributions 


-1 


cannot  be  obtained  explicitly  since  the  inverse  distribution  function  F   (•) 
cannot  be  obtained  explicitly.   An  alternative  way  of  obtaining  negatively 


correlated  (e  ,e')   begins  by  noting  that 
n  n 


e   =  K  E 
n    n  n 


and 


s  i  =  k'E' 
n    n  n 


(4.2) 


where  marginally,  from  (2.6), 
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n   n 


w.p.   (l-3)/{l-(l-a)B) 


(l-a)p   w.p.  aB/{l-(l-a)B)   , 


(4.3) 


and  E  ,  Ef  are  exponential (X)  variables,  marginally.   The  dependency  of 

n  n 

£   and  £ '  is  then  given  by 

n  n 


Cov(e  ,£»)  =  (X  2  +  C_)C__  +  [E(K)]2  C_.  , 
n  n  E  K  E 


(4.4) 


where   CT,  is  the  covariance  of  K   and  K'   and  C„   is  the  covariance 
K  n        n        E 

of  E   and  E'.   Although  any  negatively  correlated  bivariate  exponential 


n 


n 


can  be  used  for   (E,E'),  the  most  negative  correlation  is  attained  in  the 

(degenerate)  antithetic  case.   The  antithetic  choice  for   (K  ,K')>  whose 

distribution  follows  (4.1)  with  a   replaced  by   (1-B) /{ l-(l-a)B) >  does 

not  involve  degeneracy.   These  antithetic  choices  should  give  a  negatively 

correlated  mixed  exponential  pair   (e  ,e')   whose  correlation  is  almost 

as  negative  as  the  true  antithetic  bivariate  mixed  exponential  pair. 

Note  that  the  distribution  of  this  bivariate  mixed  exponential  pair 

(e  ,e')   is  a  little  complicated  in  view  of  the  break  in  form  of  the 
n  n 

antithetic  distribution  of  the  binary  pair   (K  ,K' )   at   (1-B) /{l-(l-a)B)  =  1/2 

n  n 

or  B  =  l/(l+a) .   Covariance  calculations  using  (4.4)  then  give  the  result 


_2 

X   Cov(e  ,e')  =  < 
n  n 


(l-2aB)(l-7T2/6)  -  (aB)2 


for  B  <  l/(l-kx) 


(4.5) 


(l-a)B(2-aB-B)  (1-tt2/6)-(1-B)2   for  B  >  l/(l-hx). 
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This  expression  will  now  be  used  in  determining  explicit  results  for  the 
first  autocorrelation  of  the  NEARA(l)  model.   Other  less  degenerate 
negatively  correlated  exponential  random  variables  can  be  used;  the 
simplest  and  most  easily  utilized  one  is  given  by  Gaver  (1972) . 


5.   THE  FIRST  AUTOCORRELATION 

The  first  autocorrelation  of  the  NEARA(l)  model  can  now  be  obtained, 
and  its  range  of  values  will  be  determined,  both  generally  and  in  the 
B  =  1   case,  the  so-called  TEARA(l)  model.   Interest  is  in  the  degree  to 

which  negativity  can  be  attained,  bearing  in  mind  that  with  exponential 

2 
marginal  distributions  there  is  a  theoretical  lower  bound  of  ( 1— rr  /6)  =-0.6449 

on  the  correlation.   From  (3.7),  (3.10)  and  (4.1)  we  have 


p.  =  Corr(X  ,X  .)  =  < 
1         n  n-1 


aB€-(aB) 


for  0  <  a  <  1/2 


(5.1) 


aBU-U-a)  232}/U-(2a-l)B2}  for  1/2  <  a  <  1 


This  result  is  combined  with   6   from  (4.5)  to  give  the  most  general 
expression 


P-.  = 


aB(l-2aB)(l-TT2/6)-2(aS)3, 


0  <  a  <  1/2,  B  <  l/(l-hx) 


aB  (l-a)(2-aB-B)(l-iT  /6)-aB(l-B)  -(aB)  ,    0  <  a  <  1/2,  B  >  1/U+a) 

{aB(l-2aB)(l-TT2/6)-(aB)3-a(l-a)2B3}/U-(2a-l)B  ) 

1/2  <_  a  <  1,    B  <  l/(l+a) 

{a(l-a)B2(2-aB-B)d-TT2/6)-aB(l-B)2-a(l-a)2B3}/{l-(2a-l)B  } 

1/2  <  a  <  1,    B  >  l/(l+a). 


(5.2) 
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It  is  worth  stressing  that  the  negativity  of  p   implies  the 
negativity  of  Corr(X  ,X');  then  by  virtue  of  the  general  result  (3.7) 
there  is  strong  alternation  in  the  autocorrelations  which  parallels  the 
usual   {p,)   Markov  correlation  structure  when  p    is  negative.   However, 
in  general  the  marginal  NEARA(l)  process  is  not  a  first-order  Markov 
process. 


6.   THE  LAG  ONE  JOINT  DISTRIBUTION 

Following  on  from  the  first  autocorrelation,  the  full  joint  distri- 
bution of   (X  ,X  ..  )   is  of  interest  in  describing  the  process  and  match- 
n  n-1 

ing  it  with  data.   This  joint  distribution  can  be  obtained  from  the  NEARA(l) 
model  equations  (2.8)  with  the  use  of  Laplace-Stieltjes  transforms;  thus 

*   v    (s,t)  =  E{exp(-SX  -  tX  ,)}  (6.1) 

a  ,a  _.  n     n-l 

n  n-l 

=  E(exp[-s(e   +  3V  X'  .)  -  tX    1} 
n     n  n-l      n-l 

=  E{exp(-se   -  tX  .  -  BsV  X'  _)}.  (6.2) 

n     n-l      n  n-l 


Writing     cf>      (s)      for     E(exp(-se    )}      and    taking  expectations  with   respect 

n 

to      V   ,   we  have 
n 


4>x     x        (s,t)    =   <|»e(s){a^x  x,(t,6s)   +    (l-a)<j>x(t)}.  (6.3) 

n'    n-l  ' 


The  suffix  n  has  been  dropped  from  the  right-hand  side  of  (6.3)  in  view 
of  the  stationary  assumption;  again  the  joint  distribution  of  (X,X')  is 
required.   However,  when  the  simpler  NEAR(l)  model  allowing  only  positive 

15 


dependency  is  considered  so  that  formally  X  =  X' ,  there  is  the  simpler 

result 

(s,t)    =    <j>   (s)    {a<J>x(Bs   +  t)   +   (1-a)    ^(t)}.  (6.4) 


*x,x  ; 

n     n-1 


where  4>v(t)  =  X/(X+t).   This  can  be  inverted  but  the  overall  behavior 

can  be  seen  immediately.   With  probability   1-a   there  is  a  scatter  of 

values   X  =  e    independent  of  the  X     variable,  where  as  with  prob- 
n    n  n_l 

ability  a,  X  =  e   +  BX   .   and  is  always  above  the  line  X  =  BX 

n    n     n-1  J  n     n-1 

Returning  now  to  the  NEARA(l)  model,  the  joint  distribution  of 
(X,X')   is  required.   By  constructing  Laplace-Stieltjes  transforms  from 
each  side  of  the  model  equations  (2.8),  it  follows  that 


i 

n'  n 


v  Y,(s,t)  =  E{exp[-s(e   +  0  V  X'  .  )  -  t(e  '  +  6 V'X   )  ]} 

X_ ,  X_  n     n  n-1       n     n  n-1 

=  <f>    .(s,t)  E{exp(-StV»X    -  BsV  X'   )}   .        (6.5) 
£  ,£  n  n-1      n  n-1 


Now  the  joint  distribution  of   (V  ,V)   is  available  from  (4.1)  and  so 

n     n 


4>v      y»  (s»  *-' 
n'    n 


=  4>        .(s,t) 


1   -   2a  +  acf.    (Bs)   +  a<f>Y(Bt),  0   <_  a  _<   1/2 

(6.6) 

(2a-lHx  x,       (Bt,Bs)    +    (1-a)    <j>x(Bs)    +   (1-a)    <j>    (Bt) 

n-1'    n-1 

1/2   <   a   <    1    . 


It   is   seen   that,    for     0  _<  a  _<   1/2,   <J>  ,(s,t)      is   immediately   available 

X    ,  X 
n      n 
in   terms   of  ,(s,t)      whereas    for      1/2   <^  a   <    1     a  recursive   calculation 

£  ,£ 
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is  required.   This  simplifies  somewhat  if  it  can  be  assumed  that  the  joint 

distribution  of   (e,e')   and  hence  the  joint  distribution  of   (X  ,X')   are 

n  n 

symmetric  in  s  and  t;  there  would  be  no  point  in  assuming  otherwise  for 
univariate  modelling  of  (X  }.  A  certain  amount  of  calculation  then  gives 
the  final  form  of  (6.6)  as 


[(l-2a)+a{6(Bs)  +    *    (3  t)}  ]<j>         .  (s  ,  t)  ,  0   <  a    <   1/2 

a  a  e  ,  t 


*  (s,t)    =    (  (6.7) 

(1-a)      I      (2a-l)j[4(Bj+1s)  +  *    (6j+1t)]      II      <f>         .(B^.B^), 


X    ,X' 
n     n 


j-o  x  x  i=o    £'e 


1/2  £  a    <    1. 

The  series  here  can  be  summed  in  the  3=1,  TEARA(l)  case  when  the  joint 
distribution  of   (e,e')   is  a  bivariate  exponential.   In  the  NEARA(l)  case, 
the  bivariate  distribution  which  has  been  proposed  at  (4.2)  gives 


<b    ,(s,t)  =  E{exp(-sK  E  -  tK'E'))  (6.8) 

e  ,e  n  n     n  n 


and  this  can  be  expressed  in  terms  of  the  joint  Laplace-Stieltjes  transform 

$        ,(s,t)   of  the  underlying  bivariate  exponentials.   Thus 

E ,  E 

♦e.e^8'0  =  2  1^0?  UE}E,(s,d-a)Bt)  +  ^  ((l-a)Bs, t)} 

+ 11  - 2  i=n=l)BivE'<(1-a)Bs-  (1"a)6t)  •  (6-9) 

B  <  1/d+a) 


and 
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K    B,(8,t) 


KteASit)        2   l-(l-a)B   *E,E 

1  -  2  ±-^- 

1   Z  l-(l-a)8 


[4>E  E,(s,(l-a)6t)+  <t»E  E,((l-a)s,t)], 


6  >  l/(l+a).  (6>1Q) 

The  joint  Laplace-Stieltjes  transform  of  the  distribution  of   (E,E')   in 
the  antithetic  case  is  given, with  U  a  uniform  (0,1)  random  variable,  by 

1 
$   E.(s>t)  =  E{exp[s  log  U  +  t  log  (1-U)]}  =  /   uS(l-u)t  du       (6.11) 

u=0 

which  is  a  Beta  function. 

Both  regressions  from   (X  ,X   ., )   of  the  NEARA(l)  model  are  non- 
&  n  n-1 

linear;  directly  from  the  model  equations  (2.8),  the  forward  conditional 

expectation  is 

E(X  |X  ,  =  x)  =  (l-aB)A-1  +  a3E(X*  .  |x   ,  -  x)  .       (6.12) 
n  n-1  n-I   n-1 

The  regression  on  the  right-hand  side  is  complicated  but  can  be  obtained 

in  the  6=1   TEARA(l)  case.   With  positive  dependency  only,  (6.12) 

applies  for  the  NEAR(l)  with  formally   X'    =  X   . ,  and  so  there  is  linear 

n-1    n-1 

regression  in  this  case  with 

E(Xn'Vl  =  X)  =  (i-^U"1  +  aBx.  (6.13) 


7.   THE   (X  ,X    )   JOINT  DISTRIBUTIONS;  THE  Z ^X   ,   SUM  DISTRIBUTIONS 
n  n-r  1  n-1 

These  distributions  follow  directly  from  the  basic  expressions  (3.3) 

and  (3.4)  where   X   is  expressed  in  terms  of   X     or  X'     Expecta- 

n  n-r       n-r 

tions  are  taken  over  the  independent   V  ,  V'   ,  . . .   in  turn,  with  the 

n   n-1 

following  results 

18 


<j>  (s,t)    =  E{exp(-sXTi   -    tX        )} 

X   ,a  n  n— r 

n     n-r 


r-1 

ar    n    $  (e  s) 

i-0     e 


*x  x'^'2  s^    r  odd 


{  <})X(B   s+t) 


r  even 


r-1      .  j 

+     I      aJ(l-a)      H   (f>    (BXs)    <(>Y(t). 
j=0  i=0  £  A 


(7.1) 


In  the  3=1,  TEARA(l)  case,  there  is  the  more  explicit  expression, 


*x    x      (s,t)  =  ar[*E^1-ct)s>]r 

n'    n-r 


<J>X      ,  (s,t)        r  odd 


<J>    (s  +  t)  r  even 


1  -    [a<|>   {(l-a)s}]r 

+  ci-a)  .j.  {(i-a)s>  *Y(t)1 — -r~rn — n — • 

E  X         1  -  a<f>    {  (l-a)s> 


(7.2) 


In  Section  8  these  expressions  are  used  to  derive  the  autocorrelations 
of  the  sequence  after  transformation  to  a  uniform  (0,1)  marginal  distri- 
bution. 

r 

The  distribution  of  the  sums   £ .  .  X   .   can  in  principle  be 

i=l     n-i 

obtained   from   the   expressions    (3.3)    and    (3.4)    in  a  similar  way;    for   instance, 


X     +X      ,=e     +X     ,+BVX' 
n  n-1  n  n-1  n  n-1 


(7.3) 


X+X     .+  X        =  e    +e      n+BVe'+(l  +  B2V'V)X     0  +  3V  X' 
n        n-1      n-2        n        n-1  n  n-1  n-1  n     n-2  n  n-2 


(7.4) 


Generating  functions  for  these  two  sums  can  be  written  down,  but  the  results 
get  progressively  more  complicated.   There  does  not  appear  to  be  any  simple 
general  result,  even  with  the  NEAR(l)  model. 
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8.   RUN  PROBABILITIES  AND  A  PARTIALLY  REVERSIBLE  PROCESS,  PREAR(l) 

We  have  already  indicated  in  Figures  la,  lb,  lc  that  the  sample  path 
behavior  of  NEAR(l)  processes  can  be  distinctive,  and  is  adjustable  through 
the  two  parameters  a   and  8.   This  distinctive  behavior  makes  the  model 

very  rich  and  is  principally  observed  as  runs  of  increasing  values  (up-runs) 
or  runs  of  decreasing  values  (down-runs)  or  both  (peaks) .    Such 

behavior  is  not  possible  with  Gaussian  AR(1)  models.   In  the  discussion 

which  follows  we  will  explain  the  parameterization  of  the  process 

illustrated  in  Figure  lc,  which  exhibits  a  partial  time  reversibility. 
A  simple  quantification  of  sample  path  behavior  is  given  by 

P(X  <  X  .. )  ,  which  is  related  to  the  average  length  of  up-run  sequences. 
n    n— 1 

Calculation  of  P(X  <  X  .)   follows  from  (2.3)  as 

n    n-1 


P(X  <  X  ,)  =  (1-ct)  P(X  .  >  e  )  +  aP(X  ,  >  e   +  8X  .) 
n    n-1            n-1    n        n-1    n     n-1 

=  (1-a)  P(X  ,  >  e  )  +  (1-a)  P (X    >  e  /(1-8))  . 
n-1    n  n— 1    n 


(8.1) 


By  using  the  definition  of  e    given  at  (2.6)  and  the  independence  of 

X  n   and  c    ,  the  probabilities  in  (8.1)  are  easily  calculated  and  give 
n-1        n  ° 


P(X  <  X  J  = 


1-a 


n    n-1'    l-(l-a)8 


1-6 


aS 


l+(l-a)8 


l-(l-a)6 


aB 


l+U-6)  X   l+U-8)  1(l-a)8 


(8.2) 


(1-a) (14^) 

2[l+(l-a)B]    (2-a)(l-a8) 


a(l-B) 


(8.3) 
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For  the  TEAR(l)  process  this  probability  (with  B  =  1)   reduces  to 

(l-a)/(2-a)   and  is  thus  always  less  than  one-half,  so  indicating  an 

excess  of  up-runs;  this  is  clearly  illustrated  in  Figure  lb.   A  grid  of 

values  of  this  probability  for   a, 3  =  0.0(0.1)1.0   is  given  in  Table  2. 

The  asymmetry  of  up-run  and  down-run  sequences  for  most  NEAR(l) 

processes  is  evidence  enough  of  their  irreversibility  in  time.   The  value 

of  P(X  <  X   .)   and  its  difference  from  one-half  gives  one  measure  of  this: 
n    n-1 

another  possible  measure  could  be  based  on  the  difference  between  the 

2  2 

directional  correlations   Corr(X  ,  X  .. )   and  Corr(X  ,  X  ,);  from 

n   n-1  n   n-1 

(2.3)  these  may  straightforwardly  be  obtained  as 

Corr(X  ,  X2  J  =  a3  (8.4) 

n   n-1 

Corr(X  ,  X   .)  =  a3(l  -aB  +  23)  .  (8.5) 

n   n-1 

The  equality  of  these  two  correlations  suggests  one  definition  of  partial 

reversibility,  and  for  NEAR(l)  processes  gives  the  condition  B  =  l/(2-a). 

The  simulations  in  Figure  lc  are  for  this  parametrization.   Another 

partial  characterization  of  time  reversibility  would  simply  be  that 

P(X  <  X   )  =  1/2;  surprisingly,  for  NEAR(l)  processes,  this  second 
n    n-1 

definition  also  leads  to  the  condition  B  =  l/(2-a).   Hence  we  shall 

refer  to  the  NEAR(l)  process  with  B  =  l/(2-a)   as  the  partially  reversible 

or  PREAR(l)  process.   It  is  not  fully  reversible,  even  as  far  as  the 

joint  distribution  of   (X  ,  X  n)   is  concerned,  but  it  seems  somewhat 

n   n— 1 

remarkable  that  it  is  reversible  in  both  the  run-probability  and 
directional-correlation  aspects. 
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9.   TRANSFORMATION  TO  A  MULTIPLICATIVE  PROCESS  WITH  UNIFORM  MARGINALS 

One  useful  aspect  of  exponential  processes  is  that  they  provide  a 

suitable  base  from  which  to  transform  to  other  processes  of  positive  variables; 

they  are  particularly  convenient  for  transforming  to  a  multiplicative 

uniform  process;  thus  the  transformed  process   {exp(-XX  )}   is  now  considered, 

with  derivations  of  the  autocorrelations  and  autoregressions . 

When  X   has  an  exponential  marginal  distribution  with  parameter  X, 
n 

the  variable   U  =  exp(-XX  )   has  a  uniform   (0,1)   marginal  distribution, 
n  n 

The  autocorrelations  of  the   {U  }   sequence  are  easily  obtained  from  the 

n 

joint   Laplace-Stieltjes   transform  of    the   joint  distribution  of      (X    ,    X        ) ; 


thus 


Corr(U    ,    U        )    =   {E(U  U        )    -   l/4)/(l/12) 
n       n-r  n  n-r 


12E{exp[-X(X     +  X        ) ]}    -    3 
n  n-r 


12<j)x     x        (X,X)   -    3    .  (9.1) 

n'   n-r 


Working  from  (7.1)  a  reasonably  explicit  result  for  (9.1)  is 
obtained;  the  first  expression  to  be  considered  is 


r-1 
12ar  n    <f>  (S1A)  , 
i-0  £ 


where      <}>    (s)      is    given  by    (2.5).      After   some   cancellations,    we   get 


r-1  .  r  . 

12ar     n      *    (61X)    =  6ar(l+6r)      II      {1  +   (I-oObV      .  (9.2) 

i=0     e  i=l 
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The  second  term  of  (7.1)  and  (9.1)  is 


r-1       .      j 
12(l-a)      I     aJ      n      <}>    (B1^)    <f>    (A)    -   3 
j-0         1=0  x 


r-1  .  j 

=   3(l-a)      I     aj(l+6j+1)      n    (l  +    (l-a)6i+1}"1   -   3    .  (9.3) 

j=0  1=0 


This  does  not  look  promising,  at  least  not  until  the  j  =  0   term  is  taken 
out  and  combined  with  the  -3;  the  expression  then  becomes 


r-1  j 

3(l-a)  I      aj(l  +  Bj+1)   n  (1  +  (l-a)31+1}_1  -  3a{l  +  (l-a)8}_1 
j=l  1=0 


Next   the    term     j=l      is   taken  out  and   combined  with   the   last    term;    this 
yields 


r-1  j  1 

3(l-a)      I      aj(l+ej+1)      n      {1+  (l-a)Bi+1)~1   -    3a2     n      {1  +    (l-a)3i+1}_1 
i=2  i=0  i=0 


Continuing   in   this    fashion   gives   the    final   expression 

r      r  i   -1 

-3a        n      (1  +    (l-a)3    }         . 

i=l 


Bringing  together  (9.1),  (7.1),  (9.2)  and  (9.4)  gives 
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Corr(U   ,U        )    = 
n     n-r 


6ar(14£r) 


*x  x'^x,e  ^     r  even 

l/(246r)  r  odd 


-3a ' 


n  {l+u-cOe1}  l 

i=l 


3{2(l+6rU      Y,(X,6rA)-l}      n     ( A  (r      odd) 

X'X  i=lV  l+(l-a)6     / 


=   < 


3 

2+6 


-   n     (—**-?) 

r   i=l     \l+(l-a)6X/ 


(r   even) . 


(9.5) 


This   is   the   required   result;    it   is   computationally   explicit  in  several 

cases:      the     6    =   1     TEARA(l)   model,    the  NEARA(l)    model    for  0   <  a    <^  1/2 

and   the  NEAR(l)    model    for   the   full   parameters   region.      This  latter  model 
has   as   its    transformed  autocorrelation   function 


Corr(U   ,U        )    = 
n     n-r 


2+6 


n    ( S§_^\ 

1=1  \l+(l-a)6      / 


r  =  1,2, .. . 


(9.6) 


The  only  case  of  (9.5)  which  is  not  available  in  closed  form  is  the 

NEARA(l)  model  for  1/2  _<  a  _<  1 .   The  series  expansion  from  (6.7)   for 

<$>        ,(A,6  X)   would  require  detailed  examination;  the  lower  bound  of  the 
X ,  X 

r=l   case  would  be  interesting. 

We  now  derive  the  forward  regression  E(U  |u   .. )   of  the  variables 

n  n-1 

in  this  uniform  process;  it  has  previously  been  remarked,  equation  (6.13), 
that  for  exponential  NEAR(l)  variables  this  is  linear.   As  for  the  auto- 
correlations of  transformed  exponential  processes,  equation  (9.1),  a 
a  general  result  is  available.   Without  going  into  details,  this  can  be 


written 


E[exp(-AX  )|X    =  x]*  =  A  *  A        (A,s-A)  , 
n   n-l  X  ,x 

n  n-l 


(9.7) 


where  asterisk  denotes  Laplace-Stieltjes  transform  with  respect  to   x 
of  argument   s.   Inversion  in  the  NEAR(l)  case  gives  the  desired  result 
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E(u„'Vi  =  u>  -i  1+UB)B  a-«^"e) 


(9.8) 


As  to  be  expected  it  is  non-linear.   The  corresponding  backward  regression 
is  also  available  from  (9.7). 

Finally,  we  note  that  the  results  from  Section  8  on  run  behavior 
apply  here  since  the  transformation  used  is  monotonic;  in  particular  the 
uniform  process  is  reversible  in  its  run  behavior  under  the  condition 
3  =  l/(2-a).   However,  reversibility  of  the  directional  correlations  will 
not  be  achieved  under  this  condition.   The  directional  correlations  can 
be  obtained  by  similar  methods  to  those  used  to  obtain  the  ordinary 
correlations . 


10.   ASPECTS  OF  ESTIMATION 

Formal  methods  of  estimation  are  rather  intractable  with  the  NEAR(l) 


mo 


dels:   as  an  illustration,  in  the  NEAR(l)  case  with  just  one  observation 


x  after  the  initial  value   x  ,  the  likelihood  takes  the  form 


af  (x  -3x  ),    3  <  x-j/xq  > 


L(a,B;x1,xQ)  =  (1-a)  f  (x^  + 


(10.1) 


0 


6  >  x  /xQ  , 


where   f  (.)   is  the  mixed  exponential  pdf  of  the  independent  e 
variables  given  at  (2.6).   With  more  observations,  the  full  likelihood 
becomes,  in  view  of  the  first  order  Markov  structure  of  this 
model,  the  product  of  similar  terms.   The  maximization 
needs  to  be  done  numerically  and  because  of  singularities 
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in  the  parameter  space,  the  standard  asymptotic  theory  of  maximum  likeli- 
hood is  inapplicable.   A  discussion  of  the  problems  in  the   a  =  3 
case  of  the  NEAR(l)  model  is  given  by  Raftery  (1979) .   When  a  =  1   the 
estimator  proposed  in  Gav^.r  and  Lewis  (1980)  is  the  maximum  likelihood 
estimate  (personal  communcation  from  G.  Weiss).   In  this,  section  we  limit 
ourselves  to  ad  hoc  possibilities  for  estimation  when  a  4   1. 

The  method  of  moments  can  be  developed  for  the  NEAR(l)  model:   use 
can  be  made  of  the  directional  correlations  (8.4)  and  (8.5).   The 
product  a3   in  (8.4)  is  best  estimated  by  the  first  serial  correlation, 
rather  than  the  sample  directional  coorelations .   Then  using  the  sample 
directional  correlation  based  on  (8.4)  an  estimate  of  3   can  be  obtained 
and  hence  an  estimate  of  a.   Methods  of  improving  the  efficiency  of  these 
moment  estimates  are  being  studied.   Use  of  the  run  probability  given  by 
(8.3)  is  also  a  possible  tool  for  estimation. 


11.   FURTHER  DEVELOPMENTS 

Further  work  on  this  topic  is  being  directed  at  the  estimation, 

simulation  and  sample  path  aspects.  Extensions  of  the  model  to  mixed 

exponential  variables  are  also  being  developed. 
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NEAR(l)  PROCESS— EAR(l)  CASE 
ALPHA  -  .990,  RHO  =  0.75 
BETA   -     . 758 
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FIGURE   la, 


lltli 

20  40  60  80  100 

Simulated  sample  path  for  the  EAR(l)  process  of  Gaver  and  Lewis  (1980) 
which  is  the  special  case  NEAR(l)  process  in  which  a  =  1.0.  (Simulation 
done  with  a  =  .99   to  avoid  computation  problems.)   For  this  case 

P{X   <  X   - }  =  .78   and  the  runs  of  falling  values  are  clearly 

n    n-1  °  J 

discernible. 
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NEAR(l)    PROCESS— TEAR(l)   CASE 
ALPHA    =  .758,        RHO  =  0.75 

BETA      =  .990 
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Figure   lb. 


Simulated  sample  path  for  the  TEAR(l)  process,  the  special  case 

NEAR(l)  process  in  which   8=1.   (Simulation  done  with   8  =  .99 

to  avoid  computational  problems.)   For  this  case  P{X  <X   ,}  =  0.22 

n   n-1 

and  the  predominance  of  runs  of  ascending  values  is  clearly 
discernible. 


31 


00  r 


10 


OJ 


NEAR(l)  PROCESS— PREAR(l)  CASE 
ALPHA  =  .857,  RH0  =  0.75 
BETA   =     . 875 
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Figure  lc.   Simulated  sample  path  for  the  NEAR(l)  process  which  is  partially 

time-reversible  in  that  the  directional  correlations  are  equal  and 

P(X  <  X   . )  =  1/2.   The  parametrization  for  this  PREAR(l)  process 
n    n-1 

is   8  =  l/(2-a).   Note  that  the  same  i.i.d.  exponential  sequence 
{E  }  was  used  in  the  three  simulations  of  Figures  la,  lb,  lc. 
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NEARA(l)  PROCESS-TEARA(l)  CASE 
ALPHA  =0.75 
BETA   =  0,990 
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Figure  Id.   Simulated  sample  path  for  the  TEARA(l)  process,  the  special 
case   NEARA(l)  process  in  which   3  =  1.   (Simulations  done 
with   S  =  0.99  to  avoid  computational  problems.)   Runs  of 
alternating  ascending  values  can  be  discerned,  and  are  produced 
by  the  negative  dependency  in  the  model;  this  compares  with 
the  smoother  run-up  sequences  in  the  TEAR(l)  simulation  of 

Fig.  lb. 
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